Comparison of simulated and reconstructed variations in East African hydroclimate over the last millennium

The multi-decadal to centennial hydroclimate changes in East Africa over the last millennium are studied by comparing the results of forced transient simulations by six general circulation models (GCMs) with published hydroclimate reconstructions from four lakes: Challa and Naivasha in equatorial East Africa, and Masoko and Malawi in southeastern inter-tropical Africa. All GCMs simulate fairly well the unimodal seasonal cycle of precipitation in the Masoko-Malawi region, while the bimodal seasonal cycle characterizing the Challa-Naivasha region is generally less well captured by most models. Model results and lake-based hydroclimate reconstructions display very different temporal patterns over the last millennium. Additionally, there is no common signal among the model time series, at least until 1850. This suggests that simulated hydroclimate fluctuations are mostly driven by internal variability rather than by common external forcing. After 1850, half of the models simulate a relatively clear response to forcing, but this response is different between the models. Overall, the link between precipitation and tropical sea surface temperatures (SSTs) over the pre-industrial portion of the last millennium is stronger and more robust for the Challa-Naivasha region than for theMasoko-Malawi region. At the inter-annual timescale, last-millennium Challa-Naivasha precipitation is positively (negatively) correlated with western (eastern) Indian Ocean SST, while the influence of the Pacific Ocean appears weak and unclear. Although most often not significant, the same pattern of correlations between East African rainfall and the Indian Ocean SST is still visible when using the last-millennium time series smoothed to highlight centennial variability, but only in fixed-forcing simulations. This means that, at the centennial timescale, the effect of (natural) climate forcing can mask the imprint of internal climate variability in large-scale teleconnections

Influence of fast ice on future ice shelf melting in the Totten Glacier area, East Antarctica

The Totten Glacier in East Antarctica is of major climatic interest because of the large fluctuations in its grounding line and potential vulnerability to climate change. Here, we use a series of high-resolution, regional NEMO-LIM-based (Nucleus for European Modelling of the Ocean coupled with the Louvain-la-Neuve sea ice model) experiments, which include an explicit treatment of ocean–ice shelf interactions, as well as a representation of grounded icebergs and fast ice, to investigate the changes in ocean–ice interactions in the Totten Glacier area between the recent past (1995–2014) and the end of the 21st century (2081–2100) under SSP4–4.5 climate change conditions. By the end of the 21st century, the wide areas of multiyear fast ice simulated in the recent past are replaced by small patches of first year fast ice along the coast, which decreases the total summer sea ice extent. The Antarctic Slope Current is accelerated by about 116 %, which decreases the heat exchange across the shelf and tends to reduce the ice shelf basal melt rate, but this effect is counterbalanced by the effect of the oceanic warming. As a consequence, despite the accelerated Antarctic Slope Current, the Totten ice shelf melt rate is increased by 91 % due to the intrusion of warmer water into its cavity. The representation of fast ice dampens the ice shelf melt rate increase throughout the 21st century, as the Totten ice shelf melt rate increase reaches 136 % when fast ice is not taken into account. The Moscow University ice shelf melt rate increase is even more impacted by the representation of fast ice, with a 36 % melt rate increase with fast ice, compared to a 75 % increase without a fast ice representation. This influence of the representation of fast ice in our simulations on the basal melting rate trend over the 21st century is explained by the large impact of the fast ice for present-day conditions (∼25 % difference in m yr−1), while the impact decreases significantly at the end of the 21st century (∼4 % difference in m yr−1). As a consequence, the reduction in the fast ice extent in the future induces a decrease in the fast ice effect on the ice shelf melt rate that partly compensates for the increase due to warming of the ocean. This highlights the importance of including a representation of fast ice to simulate realistic ice shelf melt rate increase in East Antarctica under warming conditions.</p

Influence of fast ice on future ice shelf melting in the Totten Glacier area, East Antarctica

The Totten Glacier in East Antarctica is of major climatic interest because of the large fluctuations in its grounding line and potential vulnerability to climate change. Here, we use a series of high-resolution, regional NEMO-LIM-based (Nucleus for European Modelling of the Ocean coupled with the Louvain-la-Neuve sea ice model) experiments, which include an explicit treatment of ocean–ice shelf interactions, as well as a representation of grounded icebergs and fast ice, to investigate the changes in ocean–ice interactions in the Totten Glacier area between the recent past (1995–2014) and the end of the 21st century (2081–2100) under SSP4–4.5 climate change conditions. By the end of the 21st century, the wide areas of multiyear fast ice simulated in the recent past are replaced by small patches of first year fast ice along the coast, which decreases the total summer sea ice extent. The Antarctic Slope Current is accelerated by about 116 %, which decreases the heat exchange across the shelf and tends to reduce the ice shelf basal melt rate, but this effect is counterbalanced by the effect of the oceanic warming. As a consequence, despite the accelerated Antarctic Slope Current, the Totten ice shelf melt rate is increased by 91 % due to the intrusion of warmer water into its cavity. The representation of fast ice dampens the ice shelf melt rate increase throughout the 21st century, as the Totten ice shelf melt rate increase reaches 136 % when fast ice is not taken into account. The Moscow University ice shelf melt rate increase is even more impacted by the representation of fast ice, with a 36 % melt rate increase with fast ice, compared to a 75 % increase without a fast ice representation. This influence of the representation of fast ice in our simulations on the basal melting rate trend over the 21st century is explained by the large impact of the fast ice for present-day conditions (∼25 % difference in m yr−1), while the impact decreases significantly at the end of the 21st century (∼4 % difference in m yr−1). As a consequence, the reduction in the fast ice extent in the future induces a decrease in the fast ice effect on the ice shelf melt rate that partly compensates for the increase due to warming of the ocean. This highlights the importance of including a representation of fast ice to simulate realistic ice shelf melt rate increase in East Antarctica under warming conditions.</p

West Antarctic Surface Climate Changes Since the Mid‐20th Century Driven by Anthropogenic Forcing

Although the West Antarctic surface climate has experienced large changes over the past decades with widespread surface warming, an overall increase in snow accumulation and a deepening of the Amundsen Sea Low, the exact role of human activities in these changes has not yet been fully investigated, which limits confidence in future projections. Here, we perform a detection and attribution analysis using instrumental and proxy-based reconstructions, and two large climate model simulation ensembles to quantify the forced response in these observed changes. We show that surface climate changes since the 1950s were driven by anthropogenic forcing, in particular the greenhouse gas forcing and stratospheric ozone depletion. Therefore, our results indicate that the 21st century changes will depend on both the greenhouse gas emissions and the ozone layer recovery

Reconciling the surface temperature–surface mass balance relationship in models and ice cores in Antarctica over the last 2 centuries

Ice cores are an important record of the past surface mass balance (SMB) of ice sheets, with SMB mitigating the ice sheets' sea level impact over the recent decades. For the Antarctic Ice Sheet (AIS), SMB is dominated by large-scale atmospheric circulation, which collects warm moist air from further north and releases it in the form of snow as widespread accumulation or focused atmospheric rivers on the continent. This suggests that the snow deposited at the surface of the AIS should record strongly coupled SMB and surface air temperature (SAT) variations. Ice cores use δ18O as a proxy for SAT as they do not record SAT directly. Here, using isotope-enabled global climate models and the RACMO2.3 regional climate model, we calculate positive SMB–SAT and SMB–δ18O annual correlations over ∼90 % of the AIS. The high spatial resolution of the RACMO2.3 model allows us to highlight a number of areas where SMB and SAT are not correlated, and we show that wind-driven processes acting locally, such as foehn and katabatic effects, can overwhelm the large-scale atmospheric contribution in SMB and SAT responsible for the positive SMB–SAT annual correlations. We focus in particular on Dronning Maud Land, East Antarctica, where the ice promontories clearly show these wind-induced effects. However, using the PAGES2k ice core compilations of SMB and δ18O of Thomas et al. (2017) and Stenni et al. (2017), we obtain a weak annual correlation, on the order of 0.1, between SMB and δ18O over the past ∼150 years. We obtain an equivalently weak annual correlation between ice core SMB and the SAT reconstruction of Nicolas and Bromwich (2014) over the past ∼50 years, although the ice core sites are not spatially co-located with the areas displaying a low SMB–SAT annual correlation in the models. To resolve the discrepancy between the measured and modeled signals, we show that averaging the ice core records in close spatial proximity increases their SMB–SAT annual correlation. This increase shows that the weak measured annual correlation partly results from random noise present in the ice core records, but the change is not large enough to match the annual correlation calculated in the models. Our results thus indicate a positive correlation between SAT and SMB in models and ice core reconstructions but with a weaker value in observations that may be due to missing processes in models or some systematic biases in ice core data that are not removed by a simple average

Role of multi-decadal variability of the winter North Atlantic oscillation on northern hemisphere climate

The North Atlantic Oscillation (NAO) plays a leading role in modulating wintertime climate over the North Atlantic and the surrounding continents of Europe and North America. Here we show that the observed evolution of the NAO displays larger multi-decadal variability than that simulated by nearly all CMIP6 models. To investigate the role of the NAO as a pacemaker of multi-decadal climate variability, we analyse simulations that are constrained to follow the observed NAO. We use a particle filter data-assimilation technique that sub-selects members that follow the observed NAO among an ensemble of simulations, as well as the El Niño Southern Oscillation and Southern Annular Mode in a global climate model, without the use of nudging terms. Since the climate model also contains external forcings, these simulations can be used to compare the simulated forced response to the effect of the three assimilated modes. Concentrating on the 28 year periods of strongest observed NAO trends, we show that NAO variability leads to large multi-decadal trends in temperature and precipitation over Northern Hemisphere land as well as in sea-ice concentration. The Atlantic subpolar gyre region is particularly strongly influenced by the NAO, with links found to both concurrent atmospheric variability and to the Atlantic Meridional Overturning Circulation (AMOC). Care thus needs to be taken to account for impacts of the NAO when using sea surface temperature in this region as a proxy for AMOC strength over decadal to multi-decadal time-scales. Our results have important implications for climate analyses of the North Atlantic region and highlight the need for further work to understand the causes of multi-decadal NAO variability

An unprecedented sea ice retreat in the Weddell Sea driving an overall decrease of the Antarctic sea-ice extent over the 20th century

Sea-ice extent is predicted to decrease in a warming climate. However, despite global warming over the past century, total Antarctic sea ice remained relatively stable from 1979 until 2015, before strongly melting. Here we explore the long-term sea ice variability by reconstructing Antarctic sea ice since 1700 CE, based on paleoclimate records and data assimilation. Our results indicate a decline in southern hemisphere sea-ice extent over the 20th century, driven by a reduction of 0.26 million km2 in the Weddell Sea that reached values at the end of the century lower than any other reconstructed period. The Ross Sea experienced an increasing sea-ice cover trend due to a low-pressure system located off the Amundsen Sea coast, offset by a decreasing trend in the Bellingshausen-Amundsen Sea. Models failed to account for the Ross Sea increase, resulting in an overly uniform estimate of Antarctic sea ice loss over the 20th century

PARASO, a circum-Antarctic fully coupled ice-sheet–ocean–sea-ice–atmosphere–land model involving f.ETISh1.7, NEMO3.6, LIM3.6, COSMO5.0 and CLM4.5

We introduce PARASO, a novel five-component fully coupled regional climate model over an Antarctic circumpolar domain covering the full Southern Ocean. The state-of-the-art models used are the fast Elementary Thermomechanical Ice Sheet model (f.ETISh) v1.7 (ice sheet), the Nucleus for European Modelling of the Ocean (NEMO) v3.6 (ocean), the Louvain-la-Neuve sea-ice model (LIM) v3.6 (sea ice), the COnsortium for Small-scale MOdeling (COSMO) model v5.0 (atmosphere) and its CLimate Mode (CLM) v4.5 (land), which are here run at a horizontal resolution close to 1/4°. One key feature of this tool resides in a novel two-way coupling interface for representing ocean–ice-sheet interactions, through explicitly resolved ice-shelf cavities. The impact of atmospheric processes on the Antarctic ice sheet is also conveyed through computed COSMO-CLM–f.ETISh surface mass exchange. In this technical paper, we briefly introduce each model's configuration and document the developments that were carried out in order to establish PARASO. The new offline-based NEMO–f.ETISh coupling interface is thoroughly described. Our developments also include a new surface tiling approach to combine open-ocean and sea-ice-covered cells within COSMO, which was required to make this model relevant in the context of coupled simulations in polar regions. We present results from a 2000–2001 coupled 2-year experiment. PARASO is numerically stable and fully operational. The 2-year simulation conducted without fine tuning of the model reproduced the main expected features, although remaining systematic biases provide perspectives for further adjustment and development.This research has been supported by the Fonds De La Recherche Scientifique – FNRS (grant no. O0100718F).Peer ReviewedArticle signat per 23 autors/es: Charles Pelletier (1), Thierry Fichefet (1), Hugues Goosse (1), Konstanze Haubner (2), Samuel Helsen (3), Pierre-Vincent Huot (1), Christoph Kittel (4), François Klein (1), Sébastien Le clec'h (5), Nicole P. M. van Lipzig (3), Sylvain Marchi (3), François Massonnet (1), Pierre Mathiot (6,7), Ehsan Moravveji (3,8), Eduardo Moreno-Chamarro (9), Pablo Ortega (9), Frank Pattyn (2), Niels Souverijns (3,10), Guillian Van Achter (1), Sam Vanden Broucke (3), Alexander Vanhulle (5), Deborah Verfaillie (1), and Lars Zipf (2) // (1) Earth and Life Institute (ELI), UCLouvain, Louvain-la-Neuve, Belgium / (2) Laboratoire de Glaciologie, Université Libre de Bruxelles, Brussels, Belgium / (3) Department of Earth and Environmental Sciences, KU Leuven, Leuven, Belgium / (4) Laboratory of Climatology, Department of Geography, SPHERES, University of Liège, Liège, Belgium / (5) Earth System Science and Departement Geografie, Vrije Universiteit Brussel, Brussels, Belgium, (6) Met Office, Exeter, United Kingdom / (7) Université Grenoble Alpes/CNRS/IRD/G-INP, IGE, Grenoble, France / (8) ICTS, KU Leuven, Leuven, Belgium / (9) Barcelona Supercomputing Center (BSC), Barcelona, Spain / (10) Environmental Modelling Unit, Flemish Institute for Technological Research (VITO), Mol, BelgiumPostprint (published version

How useful is snow accumulation in reconstructing surface air temperature in Antarctica? A study combining ice core records and climate models

Improving our knowledge of the temporal and spatial variability of the Antarctic Ice Sheet (AIS) surface mass balance (SMB) is crucial to reduce the uncertainties of past, present, and future Antarctic contributions to sea level rise. An examination of the surface air temperature–SMB relationship in model simulations demonstrates a strong link between the two. Reconstructions based on ice cores display a weaker relationship, indicating a model–data discrepancy that may be due to model biases or to the non-climatic noise present in the records. We find that, on the regional scale, the modeled relationship between surface air temperature and SMB is often stronger than between temperature and δ18O. This suggests that SMB data can be used to reconstruct past surface air temperature. Using this finding, we assimilate isotope-enabled SMB and δ18O model output with ice core observations to generate a new surface air temperature reconstruction. Although an independent evaluation of the skill is difficult because of the short observational time series, this new reconstruction outperforms the previous reconstructions for the continental-mean temperature that were based on δ18O alone. The improvement is most significant for the East Antarctic region, where the uncertainties are particularly large. Finally, using the same data assimilation method as for the surface air temperature reconstruction, we provide a spatial SMB reconstruction for the AIS over the last 2 centuries, showing large variability in SMB trends at a regional scale, with an increase (0.82 Gt yr−2) in West Antarctica over 1957–2000 and a decrease in East Antarctica during the same period (−0.13 Gt yr−2). As expected, this is consistent with the recent reconstruction used as a constraint in the data assimilation

Back to the future: Using long-term observational and paleo-proxy reconstructions to improve model projections of antarctic climate

Quantitative estimates of future Antarctic climate change are derived from numerical global climate models. Evaluation of the reliability of climate model projections involves many lines of evidence on past performance combined with knowledge of the processes that need to be represented. Routine model evaluation is mainly based on the modern observational period, which started with the establishment of a network of Antarctic weather stations in 1957/58. This period is too short to evaluate many fundamental aspects of the Antarctic and Southern Ocean climate system, such as decadal-to-century time-scale climate variability and trends. To help address this gap, we present a new evaluation of potential ways in which long-term observational and paleo-proxy reconstructions may be used, with a particular focus on improving projections. A wide range of data sources and time periods is included, ranging from ship observations of the early 20th century to ice core records spanning hundreds to hundreds of thousands of years to sediment records dating back 34 million years. We conclude that paleo-proxy records and long-term observational datasets are an underused resource in terms of strategies for improving Antarctic climate projections for the 21st century and beyond. We identify priorities and suggest next steps to addressing this.The Antarctic Climate Change in the 21st Century (AntClim21) Scientific Research Programme of the Scientific Committee on Antarctic Research are thanked for supporting the international scientific workshop at which the writing of this manuscript was initiated. This is a contribution to the PAGES 2k Network (through the CLIVASH 2k project). NJA acknowledges support by the Australian Research Council through a Future Fellowship (FT160100029) and the Centre of Excellence for Climate Extremes (CE170100023). SJP was supported under the Australian Research Council’s Special Research Initiative for the Antarctic Gateway Partnership (Project ID SR140300001). JMJ acknowledges support from the Leverhulme Trust through a Research Fellowship (RF-2018-183). FC acknowledges support from the PNRA national Italian projects PNRA16_00016, “WHISPERS” and project PNRA_00002, “ANTIPODE”. TJB, LS, and ERT were supported by the Natural Environment Research Council (NERC) as part of the British Antarctic Survey Polar Science for Planet Earth Programme. TJB additionally acknowledges support for this work as a contribution to the NERC grant NE/N01829X/1. IW thanks FAPESP 2015/50686-1, Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - Brasil (CAPES) Finance Code 001 and CNPq 300970/2018-8, CNPq INCT Criosfera 704222/200